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A topological state with simultaneous nonzero Chern number and spin Chern number is possible 
for electrons on honeycomb lattice based on band engineering by staggered electric potential and 
antiferromagnetic exchange field in presence of intrinsic spin-orbit coupling. With first principles 
calculation we confirm that the scheme can be realized by material modification in perovskite G- 
type antiferromagnetic insulators grown along [111] direction, where d electrons hop on a single 
buckled honeycomb lattice. This material is ideal for spintronics applications, since it provides a 
spin-polarized quantized edge current, robust to both nonmagnetic and magnetic defects, with the 
spin polarization tunable by inverting electric field. 
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Introduction. — The discovery of quantum Hall effect 
(QHE) by von Klitzing opened a new chapter in con- 
densed matter physics. It was revealed [Ij, 2^ that the 
integer coefficient in Hall conductance is nothing but the 
Chern number labeHng the topological property of the 
electron wavefunction. Since then to explore possible 
topological states has been one of the main driving forces 
in study of condensed matter phvsicsfsHloj. 

The breakthrough took place when Kane and Mele 
clarified that electrons on graphene, a two-dimensional 
honeycomb lattice of carbon atoms, open a gap by spin- 
orbit coupling (SOC) and achieve a topologically non- 
trivial state called quantum spin Hall effect (QSHE) Q. 
It turned out that in graphene the intrinsic SOC is very 
small which makes detection of QSHE hopeless. QSHE 
was then predicted theoretically |5|] and realized experi- 
mentally in HgTe quantum well with strong SOC [6[. 

As the common homebase for both^ Kane- Mele model 
(4*] and the spinless Haldane model [3], honeycomb lat- 
tice serves a unique role in understanding the topological 
property of electron systems. It may be illustrative to 
summarize its electronic structure paying attention to the 
Berry curvature configuration [11]. There are two sites in 
the unit cell of honeycomb structure. With nearest neigh- 
bor hopping, electronic valance and conduction bands 
touch linearly and thus form Dirac cones at two inequiv- 
alent k points, K and K', locating at the corners of Bril- 
louin zone. It is important to observe that Bloch wave- 
functions exhibit opposite chiral features around K and 
K', characterized by opposite Berry curvatures, which 
establishes the special position of honeycomb lattice in 
exploring topological state. 

With time-reversal symmetry, the Berry curvatures in 
spin-up and spin-down channels overlap as shown in fig- 
ure 1(a). Introducing intrinsic SOC opens a gap at the 
two Dirac cones, and flips Berry curvatures at K in spin- 
up channel and at K' in spin-down channel preserving 
time-reversal symmetry, which results in the Berry cur- 
vature configuration in figure 1(b), characterizing the 
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FIG. 1: (Colour online) Schematics for possible configura- 
tions of Berry curvatures and Chern numbers for insulating 
states on honeycomb lattice: (a) pristine honeycomb lattice, 
(b) QSHE, (c) trivial state with broken time-reversal sym- 
metry, (d) QAHE, (e) topological state with simultaneous 
nonzero charge and spin Chern numbers. 



QSHE state [4]. One can further flip the Berry curva- 
tures at K in both spin- up and spin-down channels, yield- 
ing the configuration figure 1(c) (a topologically trivial 
state), by a field conjugate to the valley degree of freedom 
such as the polarized field jl2l - [l5j . With an additional 
staggered electric field, which can be realized in silicene 
with buckled honeycomb lattice in terms of a uniform 
electric field[16], Ezawa could flip the Berry curvature at 
K in only one spin channel, resulting in figure 1(e) char- 
acterized by simultaneous nonzero charge Chern num- 
ber and spin Chern number [15] (see also [17l-[l9|). The 
band engineering is based on a full control on the degrees 
of freedom of spin, valley and sublattice, taking advan- 
tages of the intrinsic SOC. By the way, the Berry curva- 
ture configuration in figure 1(d) corresponds to quantum 
anomalous Hall effect (QAHE)[l8H23|. 
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In the present paper we notice that in presence of 
intrinsic SOC, a staggered magnetic field plays a simi- 
lar role as provided by the polarized light proposed by 
Ezawa[15]. Since the staggered magnetic field can be re- 
alized by antiferromagnetic (AFM) insulators, compact 
and stable devices based on the topological state in fig- 
ure 1(e) are possible as compared with the photo-assisted 
scheme. 

As material realization of our idea, we focus on d- 
electron systems in perovskite structure [24]. First, we 
choose a perovskite insulator ABO3 with G-type AFM 
order on the magnetic B atoms. As first discussed by 
Xiao et al.[9], along [111] direction B atoms form a stack- 
ing of buckled honeycomb lattice, which can be grown by 
cutting-edge molecular beam epitaxy (MBE) with atomic 
precision [25]. During the growing process, a single buck- 
led honeycomb layer of B atoms is replaced by that of 
nonmagnetic B' atoms, where the element B' is chosen 
conjugate to B in order to form a d^ configuration. For 
B'-d electrons on the single buckled honeycomb lattice, 
intrinsic SOC becomes sizable, a uniform electric field 
induces a staggered electric potential for the two sub- 
lattices, and the G-type AFM order on B atoms on the 
two sides provides an AFM exchange field. The material 
design makes the magnetic field of pure exchange char- 
acter, which avoids possible stray field from permanent 
ferromagnet. We have checked successfully our idea by 
performing first principles calculations for several mate- 
rials. In transition metal perovskites we found intrinsic 
SOC of several tens of meV, which is larger than that in 
silicene in magnitude by one order, and makes the new 
topological state available at room temperature. 

Effective model and phase diagram. — Let us illustrate 
our scheme based on the four-band Hamiltonian for elec- 
trons on a buckled honeycomb lattice under an AFM 
exchange field and a uniform electric field perpendicular 
to the plane 0] 
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XcFzTzSz + V(7z + McFzSz, (1) 



around the two inequivalent K and points where the z 
axis is taken along the [111] direction of peroskite struc- 
ture; cr's are Pauli matrices for sublattice, and = ±1 
and = ±1 are binary degrees of freedom referring to 
valley and spin. The first two terms come from the near- 
est neighbor hopping same as that in a pristine graphene 
[2^; the third term is for the intrinsic SOC with a posi- 
tive coupling coefficient A > by definition, which can be 
large due to the heavier host transition metal atom and 
the buckled structure; the fourth term is induced by a 
uniform electric potential on the buckled structure, and 
the last term is the AFM exchange field. The chemical 
potential is set to zero, and all energies are measured 
in units of the nearest neighbor hopping integral in the 
present work. As will be shown later, this model Hamilto- 
nian describes the low-energy physics of transition metal 
oxides with perovskite structure grown along [111] direc- 
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FIG. 2: (Colour online) (a) Phase diagram for staggered elec- 
tric potential V and antiferromagnetic field M, with presumed 
positive intrinsic SOC A > 0; (b) and (c) schematic diagram 
for energy level crossing upon electric field tuning for M < A 
and M > A; (d) band structure for N — 100 zigzag ribbon of 
honeycomb lattice with y = 0.25, M = 0.2 and A = 0.1 (red 
dot in (a)). 



tion. 

We reveal behaviors of the system below upon tuning 
the electric potential V while the M field and A are fixed, 
corresponding directly with our material design and ex- 
perimental manipulation. For < M < A, the gaps are 
topologically nontrivial for both spin-up and spin-down 
channels for V = 0, characterizing the QSHE state (see 
figure 1(c)). As illustrated in figure 2(b), when V in- 
creases, the gap for spin- up channel shrinks, closes and 
reopens at K' point (r2=-l) for F = A — M. This quan- 
tum phase transition turns the spin-up channel to a triv- 
ial state, while the spin-down channel remains nontriv- 
ial. This insulating state exhibiting nontrivial gap in only 
one spin channel is different from other topological states 
such as QSHE and QAHE, as discussed by Ezawa [l5|. 
When V increases further, the gap for spin-down chan- 
nel shrinks, closes and reopens at K point (r2=l) for 

V — A + M, which makes the spin-down channel trivial 
as well, and the system transforms to a trivial insulating 
state via a second quantum phase transition. 

For M > A, the system exhibits trivial gaps in both 
spin- up and -down channels for F = (see figure 1(c)). 
As seen in figure 2(c), when V increases, the gap for spin- 
down channel shrinks, closes and reopens at = — 1 for 

V = M — A, which brings the spin-down channel into a 
topological state. When V increases further, the gap in 
spin-down channel shrinks, closes and reopens at = 1 
for 1/ = A + M, which drives the system back to a trivial 



3 




FIG. 3: (Colour online) (a) A 6-Layer ABB^X perovskite with 
one layer of B (blue) replaced by B^ (grey), (b) and (c) Side 
and top views of the buckled honeycomb lattice formed by 
B^ cations, (d) Noninteracting tight-binding band structure 
of the Cg system, where each band is labeled with the corre- 
sponding electron configuration if it is occupied. Inset shows 
the two ee orbits of B^ . 



insulator state. 

Based on the above analysis we can compose the full 
phase diagram as displayed in figure 2(a) by symmetry 
consideration. The new topological state is realized when 
the absolute values of the three fields \V\ and \M\ and 
|A| form a triangle, which indicates clearly that the inter- 
vening among the three degrees of freedom, , and Sz 
is crucial for the new topological state. A typical band 
structure for zigzag ribbon of the buckled honeycomb lat- 
tice is displayed in figure 2(d); there are two edge states 
crossing fermi level which are localized at the opposite 
edges (see figure 1(e)). The present phase diagam is sim- 
ilar to the one by Ezawa [l5|, except for that the states 
with same Chern numbers lie in diagonal directions in the 
present one, while locating in same half spaces in that by 
Ezawa [l5[, due to the difference of controlling fields in 
the two Hamiltonians. 

Material realization. — In perovskite ABO3 material 
shown in figure 3(a), A and B form two penetrating sim- 
ple cubic lattices. The octahedron cage formed by six 
oxygens around the transition metal B generates the dou- 
blet Cg and triplet t2g orbits of the d electrons. In the 
[111] direction, a simple cubic structure can be viewed as 
a stacking of buckled honeycomb lattices as shown in fig- 
ures 3(b) and (c). In the new coordinate, intrinsic SOC 
appears between the two Cg orbits. In regard of material 
synthesis, recent developments in laser molecular beam 
epitaxy (MBE) technique permit one to grow perovskite 
structure along [111] direction with atomic precision (25| . 
All these make the d electrons in perovskite structure a 
very hopeful platform for realizing topological states, as 
illustrated first by Xiao et al. ^]. 
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FIG. 4: (Colour online) Band structures based on first- 
principles calculation on LaCrAgO and LaCrAuO with a su- 
per unit cell (A2B^206)/(A2B206)5 along the [111] direction. 



TABLE I: Parameters for AFM order and SOC fit from 
GGA+U+SOC calculation. For KNiInF and KNiTlF the 
electronic configurations of In^^ and TP^ are 5s ^ and 6s ^ 



field [meV] LaCrAgO LaCrAuO LaFeAgO LaFeAuO 



M 


141 


166 


541 


467 


A 


7.30 


32.91 


7.31 


33.52 


KNiPdF 


KNiPtF 


KNiInF* 


KNiTlF* 




625 


504 


290 


235 




11.38 


33.40 


5.05 


18.58 





Many transition metal oxides with perovskite struc- 
ture are known as wide-gap antiferromagnetic insulator. 
A subgroup of these materials presume the so-called G- 
type AMF order where spins align antiferromagnetically 
in the two sublattices of the buckled honeycomb lattice. 
This is the playground where we perform material de- 
sign and field manipulation to realize topological state 
with our scheme. We replace one buckled layer of B by 
a nonmagnetic element B', and obtain B'-d electrons on 
a single layer of buckled honeycomb lattice, which feel 
AFM exchange field from the G-type AMF order of the 
host material ABO3 p7j . 

There is a gap (lODq) between the Cg and t2g orbits, 
which drops t2g orbits from the low-energy physics. In 
each unit cell, two B' atoms contribute totally four Cg or- 
bits. The band structure for these four orbits is shown in 
figure 3(d) based on a tight-binding Hamiltonian with the 
hopping integrals given by Slater-Koster formula. One 
finds two flat bands and two dispersive bands crossing 
each other. Xiao et al. B [H, Mi focus on the flat 
band and developed an interesting scenario for possible 
topological state based on strong electron correlations by 
using d^(t2geg) configuration. Here we concentrate on 
the two dispersive bands by taking d^(t2geg) systems for 
which the fermi level is around their crossing point. 

In order to verify the above picture, we have performed 
first principles calculation on the host material LaCrOs 
with one layer of La2Ag206 or La2Au2 06 inserted, which 
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FIG. 5: (Colour online) Device for a "composite QHSE" in- 
sulator formed by two copies of the new topological state. 

realizes a d^(t2geg) system on the buckled honeycomb lat- 
tice. The super cell used in the calculations contains 6 
layers with a configuration of (La2B'206)/(La2B206)5. 
The calculations are performed by using density func- 
tional theory(DFT) within the generalized gradient ap- 
proximation (GGA) approach in the parametrization of 
Perdew, Burke and Ehzorhof (PBE) [28] for exchage- 
correlation as implemented in the Vienna Ab Initio Sim- 
ulation Package (VASP) [29]. The bulk lattice constant 
is set to 3. 85 A for LaCrOs and atoms are allowed to relax 
inside the super cell. The on-site Coulomb interactions 
of the 3d electrons of Cr are treated by Dudarev' method 
with an effective U value [30|, Ue//=U-J=5.4eV. We have 
checked that our main results are valid for a large window 
of the U values. The momentum space- mesh is 6 x 6 x 1 
and the cutoff energy is 440eV. 

The band structures of LaCrAgO and LaCrAuO ob- 
tained from GGA+U+SOC calculations are displayed in 
figures 4(a) and (b). First of all, the two flat bands and 
two dispersive bands close to the fermi level in figure 
3(d) are realized in the two materials. The two gaps at 
K and K', which were absent in figure 3(b) for tight- 
binding model, are induced by the AFM exchange field 
from LaCrOa and the SOC. From the two gaps at K 
and K', we figure out M = 0.141eV and A = 7.3meV 
for LaCrAgO, and M = 0.166eV and A = 32.91meV for 
LaCrAuO. We have also calculated several other possi- 
ble materials and the results are summarized in TABLE 
1. As the maximal nontrivial gap is given by 2 A, the 
results from first principles calculations suggest that a 
topological material may be realized at room tempera- 
ture. Typical electric field to realize the new topological 
state in the order of O.IV/A. 

Possible applications. — This new topological state 
available probably at room temperature is ideal for spin- 
tronics appHcations, since it provides a spin-polarized 
quantized edge current, robust to both nonmagnetic and 
magnetic defects, with the spin polarization tunable by 



inverting electric field. 

Using two copies of the new topological state, one can 
realize a "composite QSHE insulator". As shown in fig- 
ure 5(a), we fabricate two pairs of electrodes (golden) 
with opposite bias voltages. The two patches exhibit 
states with opposite chiralities and spins. Spin Hall con- 
ductance can be measured based on this device. As 
shown in figure 5(b), when a positive spin current Is 
is injected from the left electrode, the out-going spin- 
up current flows from the left electrode along the top 
edge of upper patch (upturned red solid curve), while 
the in-going spin-down current flows into the left elec- 
trode along the straight line between the two electrodes 
(solid green line), which results in zero charge current. 
Because charge current flows only on the top edge and 
cumulates charges there, a Hall voltage drop can be de- 
tected between the top and bottom edges of the de- 
vice. As compared with the QSHE realized before, even 
a time-reversal-symmetry-broken perturbation, e.g., a 
magnetic field or magnetic impurities cannot backscat- 
ter the counter charge currents with spin up and down 
at the center of the device, since the charge currents are 
now distributed in different patches as shown in figure 
5(b). 

Discussions. — The present topological state is robust 
against Rashba SOC [H, [23|, which, although small, 
should be present since the electric field breaks the in- 
version symmetry with respect to the honeycomb lattice 
plane. Actually, it is known that in graphene the non- 
trivial band gap closes for Xr > X and the QSHE will be 
destroyed [4] due to the mixing of the two spin channels by 
Rashba SOC. In the new topological state with for exam- 
ple F = M > 0, a gap of 2(2M — A) is opened in spin- up 
channel whereas the gap in the spin-down channel is 2 A. 
Since M ^ A as shown in Table 1, even a Rashba SOC 
stronger than the intrinsic one Xr ^ X cannot close the 
gap by mixing the two spin channels. 

Note. — About finishing our manuscript, we became 
aware of a recent work by Ezawa arXiv. 1301. 0971, ad- 
dressing the same physics. In his work, the staggered 
magnetic field is generated by two ferromagnets on the 
two sides of silicene layer. 
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